Metastable states of a driven flux lattice in a superconductor with strong pins 
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The flux lattice driven by a uniform driving force in a superconductor with hot, strong, sharp 
and randomly distributed pinning centers, with applied magnetic field half the matching field is 
simulated. At low temperature both a non activated regime, where flux motion occurs within a 
robust percolative flux flow channel, and an activated regime are obtained depending on the sample 
preparation. These two regimes exhibit distinct resistivity and magnetic induction. In the non 
activated regime, a clear flngerprint is observed in the autocorrelation function of the longitudinal 
resitivity, which oscillates at a frequency close to the inverse lattice diffusion time. 
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I. INTRODUCTION 



High-temperature superconductors have revealed a 
complex static vortex lattice phase diagram in the 
temperature-magnetic field (or disorder strength) plane, 
and it is now well established, both theoretically and 
experimentally, that at least three distinct phases are 
present : the high temperature liquid phase, the quasi- 
ordered and the glassy solid phases at low tepjiperature 
for low and high magnetic field respectivelylJTj. Intro- 
ducing a driving force as a third axis, leads to the even 
more complex dynamic phase diagram of a driven vor- 
tex lattice. A theoretical insight on this situation has 
been given in refs. ^,||,^-10. 



The emerging picture is 
the existence, for innceasing drive, of creep, plastic and 
moving solid regimeaS. The moving solid, once thpught 
a 'moving crystal' within the perturbation theoryQ, was 
found to be a 'moving glass^p, which vortices flow along 
elastically coupled channelsBll3. Such a theoretical work 
has widely benefited from the comparison with numer- 
ical experiments. Simulations first demonstrated the 
existence of a plastic flow xegupe within channels for 
the highly defectija^atticeO'EJ'EJ, yielding non-linear 
I-V characteristicsd'Ej. Dynamic phase transitions or 
crossovers between the plastic regime and moving solids 
for lattices driven by a uniform external force were put 
into evidence in refs [7|Jl^-^. The case of the field gra- 
dient driven lattice (Bean state) was extensively studied 
by F. Nori et al (see e.g. ref. ^ and refs therein). Be- 
sides these effort to determine a univoque dynamic phase 
diagram for the driven lattice, some evidence exist for 
the existence of multiple metastable states. Indeed, the 
static glass is characterized by the existence of infinite 
barriers that prevent the exploration of the entire phase 
diagram and the existence of a degenerated ground state. 
As a consequence, one might also expect that the driven 
lattice can show multiple metastable states close to the 
static glassy phase, i.e. in the plastic regime. Such a 
behavior was observed in ref. where it was shown at 
zero temperature that some filamentary flow channels are 
stable in a finite range of driving force. The transition 



between different flow channels structures as the driving 
force is varied was shown to result in steps in the I — V 
curves. In the following, I show that a random distribu- 
tion of strong pinning centers allows for the realization of 
metastable states of the disordered driven lattice at non 
zero temperature. 



II. EXPERIMENTAL DETAILS 

A two dimensional lattice submitted to a uniform driv- 
ing force in the presence of 'hot pins' is simulated. A 
semi-infinite sample is considered. Two opposite edges 
of the sample are submitted to the external magnetic 
field, -Bq, which is simulated by an extra force /s^ acting 
on each vortex, perpendicular to the sides of the sample. 
The force acting on a vortex at a distance x from the edge 
is the one imposed by a semi-infinite vortex lattice at a 

1/2 

distance -I- x, where ao = {^a/ B) ' is the flux lattice 
spacing at the equilibrium. Whenever the force acting on 
a vortex situated at the sample edge is directed towards 
the inside of the sample, this new vortex is introduced. 
The other two edges of the sample are subject to a peri- 
odic boundary condition, thus simulating a semi infinite 
geometry. 

Flux lines are assumed rigid rods. This correctly mod- 
elizes a layered superconductor with decoupled layers, or 
a three dimensional material with a large line energy. 
The force per rUliit length between vortices separated by 
a distance r isc3: 



(1) 



where K\ is a Bessel function. This is a good approxi- 
mation strictly only in the case of vortex lines (rods) and 
for 2D vortices a logarithmic interaction should be used. 
In the present case, the more rapid decrease of the Bessel 
function allows us to cut the interaction between vortices 
at a distance 5A and save computation time. 

Strong pinning centers are randomly distributed in the 
sample. The density of the pinning sites is n = -B^/^&Oj 
with $0 the flux quantum and the 'matching field' 
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for which an equilibrium flux line lattice shows the den- 
sity of flux lines n. The pinning sites are assumed normal 
cylinders parallel to the applied field, with radius cq. The 
situation where, at low temperature, the vortex core ra- 
dius ^(T) is smaller than cq and pinning is due essentially 
to the reduction of the core energy when the line sits on 
the pin is considered. The force per unit length exerted 
by a pin at a distance r to the line is given by: 

/p(r) = aeorjrQ^ for r < for r > (2) 

where £o = ('&o/47rA)^ is the line energy, rg = cq -I- 'f/2 
and (T < 1. The pinning force is exactly balanced by the 
Lorentz force, J A #o, for J = Jp = ctEo/C'^'o- From 
eq. |l| and |^, the energy scale for vortex repulsion is 
\Jvv — '2£oKo{ao/X) and the one for pins attraction is 
Up = creoro/2^. 

The temperature in the whole sample is assumed zero, 
excepted on the pinning sites ('hot pins' model). That is, 
thermally activated depinning from the pinning centers is 
present, while thermal fluctuations of the unpinned vor- 
tices are absent. Within such a model, there is no melting 
of the flux line lattice, and the contribution of the ther- 
mal fluctuations of the unpinned vortices to depinning is 
neglected. During the simulation, the force acting on the 
vortices, and the corresponding activation barrier for the 
lines which are pinned, are computed. The equation of 
motion: 

r] dr/dt = fBo + fvv + fp + J A 00 (3) 

is applied for a time step At small enough so that vortex 
displacements within this interval are much smaller than 
Oq. For the pinned vortices, escape from the wells pro- 
ceeds according to the Arrhenius exponential probability. 
During the simulation, the vortices trajectories are reg- 
istered, as well as the resistivity, normalized to its flux 
flow value : 

p/pF - (JBoAt)-^ ^ + At) - x,{t) (4) 

i 

In the rest, the driving current is normalized to Jp and 
the thermal energy fe^T to the pinning energy Up. A 
sample of dimensions 40ao along the driving current and 
SOflQ perpendicular to it was studied. The following 
superconducting parameters were used: A = 1400 A, 
^ = 18 1, ro = 2^, cr = 0.1, yielding Up/Uy^ = 0.075, and 
B(j, = 2Bo — 5000 G. This situation, using the terminol- 
ogy in ref. ^2], is the one of random, sharp (rp ^ ao) and 
strong {Up/roCee ^ 1) pins close to the matching field. 

III. RESULTS AND DISCUSSION 

I have first computed the resistivity and the induc- 
tion at T = using two different procedures (Figj^). 
Both procedures start from a regular hexagonal lattice at 
equilibrium with the external magnetic field and J — 0. 



In the first procedure, the driving current is first set to 
J — Jp. It is then decreased by steps down to J = 0. 
A steady state - characterized by a steady resistivity 
and induction - is allowed to settle at each step. In 
the second procedure, the driving current is directly set 
from zero the measuring value. I denote A and B these 
two procedures, and pA and pB the associated resistiv- 
ity. From these experiments, it is possible to distinguish 
three regimes. I: for J < 0.1, pA = and pB — 0. II: 
for 0.1 ^ J ^ 0.5, PA and pB is found randomly 
on two branches, one of them coinciding with pA, the 
other being zero or close to zero. Ill: for 0.5 ^ J, both 
procedures yield similar non zero values. The existence 
of a clear transition between regimes II and III is also 
clearly put into evidence by procedure A only, as both pA 
and the induction show a marked discontinuity (FigJ^). 
In regime II, the existence of two branches for proce- 
dure B mirrors in the induction which also splits into 
two distinct curves, one of them coinciding with the one 
obtained with procedure A. For all regimes, no Bragg 
peaks could be detected from the time-average structure 
factor and flux lines ordering is amorphous. Vortices tra- 
jectories in regimes II and III both conai£i_i|n irregular 
channels, characteristic of plastic motiontil^t^. 

I now argue that the upper branch of the resistivity in 
regime II is anomalous with respect to thermal activation 
for depinning. Heating a sample with 0.1 < J < 0.5 on 
the p = branch from T = to finite values yields an ex- 
ponentially activated resistivity (Fig. The activation 
energy is found up to one order of magnitude lower than 
the one in a one-dimensional model, U{J) = (1 — J)^, 
which is due to the inhomogeneous screening current in 
two dimensions. In a striking different way, the upper 
branch of the resistivity is found stable when tempera- 
ture is increased, until the system switches to the acti- 
vated branch. Comparison of the vortices trajectories for 
the activated regime and the non-activated one shows lit- 
tle difference (Fig. H), but the autocorrelation function 
of the resistivity noise, < p{t)p{t + t) > t — < p >^ , 
exhibits oscillations at a frequency close to the inverse 
lattice diffusion time, 1/rao — <hiJ /c rj ag, in the case of 
the non-activated regime (Fig.g). This is the signature 
of the coherent motion of vortices within flux flow chan- 
nels, whereas activated motion tends to randomize flux 
line motion. The distribution of the activation energy 
for single vortex depinning, as shown in Figjo] (we are 
here only concerned by individual depinning, inherent to 
our thermal model), provides further indication that the 
upper branch of the resistivity could be stable against 
thermal depinning. It is clear from this data that, as 
one enters regime II, the distribution splits into a dirac 
function at U — for the unpinned vortices, and a wide 
component that peaks at a non zero value with a gap at 
low energy. It is then appealing to credit the observed 
stability at low temperature of the pinned 'edifice' to the 
existence of this gap. On the other hand, I have observed 
that two different runs where the sample is heated from 
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the p 7^ branch can yield two shghtly different tem- 
peratures for the switch to the activated regime. Also, 
as can be seen in Fig. ^, the system may spontaneously 
switch from the apparently stable non-activated regime 
to the activated one at a constant temperature. Whether 
the flux flow channels are truly stable upon thermal ex- 
citations or there is an increase beyond the experimental 
window of some relaxation time is thus an open question. 

Clearly, the existence of regimes II and III is related 
to the possibility for the flux lattice to store elastic en- 
ergy (proportional to B — Bo in Fig^ in the compression 
(II, B > Bq) or the extension mode (III, B < Bq), as 
the transition between these two regimes coincides with 
the one between these modes (FigJ^) . In the same way, in 
regime II, the unactivated branch involves the compres- 
sion mode, whereas the activated one involves the exten- 
sion mode, as in regime III. This is also clearly seen in 
Fig.^, where bursts in the resistivity after the system has 
switched to the activated mode correlate with the ones 
of the induction. Then, both the existence of two dis- 
tinct regimes and of the anomalous unactivated branch 
for the resistivity are a consequence of the applied field 
being close to the matching field: for lower ratio B/B^, 
I expect regime II to shrink. 

Finally, the characteristic frequency for correlated mo- 
tion within flux flow channels may be evaluated from the 
conventional Bardeen-Stephen expression for flux flow re- 
sistivity, yielding I/tqo = Jp„c/aoi/c2- I consider the 
case of 277-NbSe2 which is a good candidate for the ex- 
perimental realization of the above situation, as it ex- 
ibits low depinning critical current for as-grown crystals, 
as well as a small line energy. We first need to evaluate 
the depinning current at the transition between regimes 
I and II, Ji/n- At this point, the driving force on the 
flux flow channel is balanced by the shear force exerted 
on it by the pinned lattice : Jj 



tocorrelation function of the longitudinal resistivity at a 
frequency close to the inverse lattice diffusion time. 
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eo/4$o cLo- Using parameters in section^ 
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we 



obtain J//// — 0.05 Jp, in reasonable agreement with 
the experimental data (Fig. |l|). Going back to the case 
of 2^r-NbSe2, using Hc2 = 2 10^ Oe, A = 2000 
Pn = 5 10~^ cm for the normal state resistivitytfta, 
one obtains Jj/n « 2 10^ B^/^ A cm"^ ^nd I/tqo ^ 90 B 
kHz. For induction larger than a few lO'^ G, Jj/n is 
well above the critical, current density usually observed 
for as-grown crystals^, so that pinning by natural de- 
fects should not prevent the observation of flux motion 
within the channels. However, the characteristic fre- 
quency strongly exceeds the bandwidth of conventional 
four probes transport measurements. 

In conclusion, it was shown that a random distribu- 
tion of pinning centers with applied field close to the 
matching field allows for the preparation of two compet- 
itive regimes of the disordered driven fiux lattice, within 
a large range of driving force. One of them involves 
non activated fiux fiow motion within channels, which 
are found robust against moderated thermal depinning. 
This regime shows characteristic oscillations of the au- 
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FIG. 1. Resistivity and magnetic induction for procedures 
A and B as described in the text {T — 0). 
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FIG. 2. Vortices trajectories (J = 0.2) at ksT = 0.04 (a) 
and fcsT = 0.025 (b, c). The resistivity is the same for con- 
figurations a and b. During the recording, 90 vortices have 
crossed the sample edge for all three samples. For sample b, 
motion within the channels is non-activated. 
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FIG. 3. Arrhenius plot of the resistivity upon heating from 
the p — branch in regime II in Fig.^ (dashed line) and from 
the p 7^ one (full line). The dashed line may be described 
for both increasing and decreasing temperature. 
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FIG. 4. Resistivity and induction at fcsT — 0.03 showing 
a spontaneous switch from the unactivated regime to the ac- 
tivated one. The sample has been progressively heated from 
the p / branch in Figjl] (J = 0.2). 
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FIG. 5. Logarithm of the autocorrelation of the resistiv- 
ity noise for the three samples in Fig.^. The autocorrelation 
function for the unactivated sample shows marked oscillations 
at a frequency « 0.7 the inverse lattice diffusion time. 
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FIG. 6. Distribution of the activation energy, i7, at T = 
(normalized to area unity). Below J ~ 0.5, the distribution 
for the pinned vortices shows a gap at low energy. 
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